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I . INTRODUCTION 


The solar collector discussed In this paper is an all-glass, selectively 
coated, evacuated collector. The collector consists of three glass tubes 
inside one another as shown schematically in Fig. 1. The middle tube 
is the absorber tube which is covered on the outside with a coating that 
has the properties of high absorptance in the solar spectrum and low 
emittance in the infrared spectrum. The largest tube or cover tube surrounds 
the absorber tube and is sealed to it at one end. The annular space between 
these two tubes is evacuated so that heat transfer by conduction and con- 
vection is essentially eliminated between the two tubes. At the same time 
the cover tube provides protection from the environment for the selectively 
coated surface. A third, smaller tube is inserted inside the absorber tube 
for delivery of the working fluid; thus the inlet and the outlet for the work- 
ing fluid are at the same end of the tube. This design makes manifolding of 
the collector tubes into an array convenient. Table 1 gives the dimensions 
of the three tubes in a single assembly, while Fig. 2 is a photograph of a 
collector array at the Owens-Illinois, Inc. test site in Toledo, Ohio. 

II. COLLECTOR PERFORMANCE MODEL 

The useful heat obtained from a single collector tube can be expressed 
in a form similar to that for flat plate collectors [1]: 

«„ = (’) 
where A^ is the absorber tube diameter times the collector length and Aj^ 
is equal to ttA . 

Three terms in this equation are discussed in some detail in the 
following paragraphs: Fj^, the performance index; the effective 

insolation on the collector; and U^, the collector loss coefficient. 
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A. Collect cf-r Performance Equation 


Equaticsjs (2) and { 3 ) are derived 


from the heat balance of control volumes involving the delivery tube, 
annulus and associated tube surfaces. Figure 3 indicates the location 
of these control volumes. Since water is assumed as the working fluid, the 
heat capacity of the several glass tubes can be neglected as small compared 
with the heat capacity of the water. 

PfCpA^i(9Tj/8t)+mCp(3T^/ax)+UjPjT^-UjP = 0 (2) 


Pf'^p\2 


3T^/9t)-mCp{8y3x)-U^PJj+{UjP^+ 


UlPl/[1 + (UlPl/^ 3^3^^^'^2 = {l/[l+{ULPL/U3P3)]}(“TSg^fPc+ULPLTa) (3) 


The equations (2) and (3) were solved in terms of the following reduced 
variables : 

5 = (UjPj/BCpjx ; t = (UjPj/PfCpA^^)* 
and the following parameters: 


» = i b' = b * 1-F' 

F' = 1/[1 + (U^^P|_/U3P3)] ; * ’'a 

A as 

The steady state condition occurs when aTj/st = 3T2/3t = 0 and T^ is 
constant. These conditions never exist in practice, but in many situations 
the collector can be described as operating in a "quasi" steady condition 
where the steady state solution gives a good approximation to the collector 
behavior. 

With the steady state assumptions, and the boundary conditions: 

Ti(0) = T^^ and = T 2 (Ci), (5i = UaPat/mCp), the fluid temperatures 

along the length of the collector tube are given by equations (4) and (5). 

Ti(0 = y(T.^-Tg)e“i^ {Ccoshaj2{5i-5)+(u),i/u)2)sinhw2(Ci-5)]/ 

[C0Shu)2? l"*" (bu l/^2 1 3 1 
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1-5 )- (0)2/102 )sinhu2 (Ci~C) 3/ 

[cosho)25i+(oji/a)2)sinho)25i]} (5) 

vvhere 

0)1 = UlPlF72U3P3 

0)2 = UlPlF’[1+4UiPi/Uj_PlF’]^/72U3P3 
The useful heat obtained from the collector tube is: 

Q^i = mCp[T2(0)-Ti(0)] = {sinho)2Ci/[o)2?i(cosho)2ei + 

(o)i/o)2)sinho)2^i)]}F'A^[aTSg^^-U^_(A|_/A^){T^.^-T^)] f5} 

Equation (7) provides the definition of Fj^/F', the performance index. 

Fj^/F' e sinho)2Ci/{o)2?i[cosho)2Ci+(ioi/o)2)sinhoj2?i]} {'?) 

B. Loss Coefficient The effectiveness of a solar collector is 
determined primarily by the amount of heat lost from the collector during 
operation. This heat loss governs the potential operating temperature of 
the collector and the level of insolation required in order for opera- 
tion to be feasible. Fundamentally, the heat loss is the product of a loss 
coefficient, loss area and temperature difference. While the heat loss can 
be reduced by reducing any one of these terms, the loss coefficient is the 
term that is most readily made smaller by application of engineering and 
scientific techniques. Various schemes have been proposed and used to re- 
duce the loss coefficient In this collector, a vacuum and spectrally 
selective coating have been used to good effect. 

The heat loss in the collector is controlled by radiation loss from 
the selectively coated absorber surface to the cover tube inner surface. 
Using the subscript scheme in Figure 3 to identify the surfaces, a com- 
bined radiation, conduction and convection heat loss network is drawn 
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in Figure 4. 

The heat lost by the collector is given by Eq. ( 8 )» 

Ql = ULA4(T4-Tg) (8) 

v/here 

Ul = [1/hi + l/h2 + l/ha]-! 
and 

hi = cr(T4+T5)(T42+T52)/{[(l-E4)/^4]+(l/F45)+[(l-e5)/e5](A4/A5)]} 

h 2 = k/[(D 4 / 2 )ln(D 6 /D 5 )] 

h3 = [h+E6a(T6+Tg)(T62+Tg2)](A6/A4) 

The loss coefficient is that for an area A 4 , which is the col- 
lector loss area A|^. 

C. Insolation on Tubes in an Array The collector tubes are axially 
symmetric and iiave an aperture of 360° for light collection, facts that can 
be used to advantage in building up multi-tube collector arrays. As will 
be shov/n in this section, the insolation available to a tube can be increased 
substantially if the tubes in an array are suitably spaced apart and a re- 
flecting screen is placed behind them. As a result, it is possible to 
achieve a highly cost effective collector array whose thermal performance 
on an installed area basis is excellent. 

Spacing of the tubes enhances the available insolation in a number of 
ways, First, each tube has an intercept area which can be made independent 
of sun angle between limits set by tube spacing, v/hich determines when 
shading from neighboring tubes occurs. Second, the diffuse light available 
to a given tube is generally greater than that available to a planar sur- 
face occupying the same project area, and this component increases with 
tube spacing. Finally, light passing through the gaps between the tubes is 
reflected onto the undersides by the backing screen, and the amount intercepted 
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by a tube Increases with tube spacing. 


ORIGINAL PAGE IS 
OP POOR QUALITY 


The details of these effects are discussed in the follovring paragraphs 
for the case In which the tube axes are In a north-south orientation. 

While either diffusely or specularly reflecting surfaces In a variety of 
shapes can be used as a backing screen, the case to be examined Is that 
of a planar, diffusely reflecting surface such as Is used with present 
tubular collector arrays. A major argument for choosing such a surface 
Is that It could consist of standard roofing materials over which the 
array would be Installed, and thereby enhance performance significantly 
while adding little or no cost to the collector. All that would be 
needed Is a reflective outdoors paint to Increase the reflectance of the 
screen. , 

1. Beam Component A tube In an array Intercepts the beam component 

I 

of Insolation both directly, on the cross sectional area of the absorber 
tube, and Indirectly from reflections off the backing surface. The directly 
Intercepted beam component can be determined by considering a tube axis In 
a north-south orientation tilted at an angle s above the horizontal at 
latitude L. If Sgg is the beam radiation in a plane perpendicular to the 
sun's rays, which are incident at declination 6, and the hour angle u Is 
taken to be zero at solar noon. It is straightforward to show that the 
component Sgg Intercepted directly by the absorber tube is 

S'bd = SgQ{l-[sin(s-L)cos6cosa)+cos(s-L)s1n6]2}i/2 (9) 

Equation (9) Is valid only so long as the tube is not shaded by neighbor- 
ing tubes. For an array with tubes spaced at distance d, however (see 
Fig. 5), shading begins to occur at hour angles |to| > |o)q|, where 

|o)ol = cos- 1[(D4 + D6)/2d] (10) 
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Shading can be taken into account through a factor g{toj, wnich is given by 


g(w) = 1 , |ul s |ojq| 

d } I Ds\ 

~ COSu + •- I 1 |»|‘^l^l<^nl 

D4 2 \ D4/ 

The beam component intercepted directly by a tube in an array is then 
SgD = quantity Sgp can be written in terms of the beam 

radiation on a horizontal surface, Sgj^, by using the relation Sgg = RjSg^, 
where [2] 

P _ / V {l-[sin(s-L)cosi5cosa)+cos(s-L)sin6]^}^/^ 

cos6cosLcosd)+sin6sinL (12) 

In practice, it is a good approximation to assume that D4 Dg for 
the purpose of evaluating the shading factor g(aj). Then, for example, if 
d = 2D0, |o)q| = 60° and the individual tubes are unshaded for eight hours 
a day. If the tilt angle s is equal to the latitude L, it can be seen 
from Eq. (9) that Sgg is independent of oj during this eight hour period. 

In fact, one finds in this case that 

SgQ - SgQ cos6 , 1 0) I ^ cos“^Dg/d 

SgQd 

= COS6COS0), U| > cos"^Dc/d 

Dg (13) 

In addition to the directly intercepted beam component Sgg, the tubes 
in an array also receive a component Sg^ due to reflections of beam light 
from the backing surface. In general, the back-reflected light will have 
specular as well as diffuse character, but for the purpose of estimating 
the component Sg^ we will assume completely diffuse reflections. Figure 
5 specifies the geometry of the array. Beam radiation passes between the 
spaced tubes, giving rise to a series of light strips whose widths W and 
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center positions relative to the axis of a given tube both vary with 
hour angle w. 

W = d - De/coso) (14) 

= Dgtano) + (i + l/2)d (15) 


The rate at which beam energy arrives at each strip is = RpSg^W' 2 ., 
where Rp is the angle factor for converting beam radiation on a horizontal 
surface to a south facing surface tilted at angle s [2]: 

P _ cos (L-s)cos6cosa)-<-sin(L-s)sin6 

T ~ cosLcos6cosw+sinLsin5 (16) 

If the backing screen has diffuse reflectance p, the rate at which energy 
is reflected is = pRpSgp|Wji. The fraction of the energy reflected 

from a strip at X^.j that is incident on the tube of interest is therefore 


Qgy(i) pRpSg,^W£F.^ 

where F^j is the geometric shape factor of the strip at X^.j and the tube. 

The total energy Qgj from all strips can be obtained by summing the contri- 

butions in Eq. (17) over all strips. In addition to the light from the 
strips, a given tube will also receive reflected light from the overhanging 

screen at either end of the array. While this contribution can be important 

for tubes near the ends of the array, it is neglected here both for simplicity 
and because its importance to the overall performance of a large array is 
small. Thus, if the back-reflected beam flux is defined as Sgp = Qgj/D^^ 
(i.e., energy per unit of absorber tube cross section area, for consistency 
with the directly intercepted flux, , one has 


’BR 


W_ 

D4 


pRrSbh 


•"iT 


(18) 
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The shape factors for those strips whose views of the tube is 
unobstructed by other tubes can be evaluated by Hottel's crossed and 
uncrossed string method [3]. In fact, it is found that 


Dit W/Dr 

F-t = — tan-i 5 

^ ' 2 W 1 + ( X .^2 . w 2 / 4 )/ d ^2 


( 19 ) 


For those strips whose view is partially obstructed by other tubes, the 
shape factors are more complicated. However, such strips lie relatively 
far away from the tube of interest so that the contributions of their shape 
factors to the sum in Eq. (18) is relatively small. Thus, it is a good 
approximation to use Eq. (19) for all the shape factors appearing in Eq. 
(18). Combining Eqs. (18) and (19), becomes 


’BR 


!!e!m j tan-i 


W/D 


B 


2 i 1 + (X..r^ - WV4)/Dg2 (20) 

Equation (20) is cumbersome to evaluate and does not lend itself 
to simplified approximations for arbitrary tube spacings d and screen 
distances Dg. A detailed analysis shows that for d 5 2De an excellent 
approximation to Eq. (20) is 


where 


^BR ” 
De 


( 21 ) 


A H 


D4 i 


, ^-T 


o )=0 


( 22 ) 


For spacings d > 2De, however, Eq. (21) is not a good approximation to 
Eq. (20). 


The portion of the beam component striking the outer edges of the 
cover tubes is near grazing incidence and therefore is largely reflected. 



Some of this radiation reaches neighboring tubes both directly ana in- 
directly, via secondary reflections from the screen. Test results have 
proven the overall contribution of these reflections to be small, and 
so it is neglected in the present analysis. 

2. Diffuse Component As in the case of the beam component, the 
diffuse component of insolation is intercepted both directly by the tubes and 
indirectly from reflections off the backing screen. If the apparent origin 
of the diffuse radiation is localized over a region of sky near the solar 
disk, as might be the case on a clear day, it is reasonable to treat the 
diffuse component as beam radiation and use the total insolation in place 
of the beam insolation in the equations developed in Section C.l. If, at 
the other extreme, the diffuse component is distributed uniformly over 
the sky dome, as might be the case on a cloudy or hazy day, its contribu- 
tion to the total insolation on an array of tubes must be calculated 
separately. In the discussion that follows, it will be assumed that 
the diffuse insolation is distributed uniformly over the sky dome. 

Assume first that the collector tilt s is zero so that ground re- 
flections need not be considered. Then, if is the diffuse insolation 
from the sky dome and the shape factor of an absorber tube and 

the sky, the diffuse flux intercepted directly by the absorber tube is 

^dD " ’^'^TS^d (23) 

As before, the flux ^dD defined per unit of cross section absorber 
tube area. 

The component of diffuse radiation reflected off the backing 
screen must be added to SdD- Consider a thin strip of backing screen with 
width dx and long dimension a parallel to the tube axes. Let ^(G) 
be the shape factor of this strip to the sky through a particular gap G. 

Then the total view factor of the strip to the sky through all gaps is 
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2 

^dx.S " G (24) 

where the sum includes all gaps G through which diffuse light can reach . 
the strip dx. The back reflected radiation per unit of absorber tube 
cross section area is therefore 


dSdR = ^pS/dx.S^’^T.dx (25) 

where dFj is the shape factor of a given absorber tube to the strip 
dx and p the reflectance of the screen for diffuse radiation, here assumed 
equal to that for the beam radiation. Integrating over all strips dx, 

^dR ^ ^P^d •Ccreen^''T,dx’'dx,S ~ (26) 

where Fjp is the shape factor from an absorber tube to the backing screen 
and the function F is defined as 


^ ~ Fjp *'screen‘^^T,dx*^dx,S 

For a large enough screen, it is a good approximation to take Fjp 
so that the total diffuse insolation on an absorber tube becomes 



SdT = ^ 

The factor Fj^ is readily evaluated by Hottel’s crossed and uncrossed 
string method [3] and is found to vary with tube spacing d between the 
limits 0.27 5 Fjg ^ 0.5. The lower limit corresponds to close packed 
tubes (d = De), while the upper limit corresponds to infinitely spaced tubes. 
The function F defined by Eq. (27) depends on the shape factors F^^ ^(G) 
through Eq. (24). Each F^^ ^(G) gives the intensity distribution of 
light on the screen resulting from diffuse radiation passing through 
a gap G, and by using the crossed and uncrossed string method [3] it 
can be shown that in general. 
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F+(G) 


(d/2 ± x)[(d/2 ± x)2 + Dg2 - D62/4]V2 - DgDe/2 

(d/2 ± x)2 + Dg2 


( 29 ) 


where the origin of x is at the center of the gap. 

When the collector tilt s is other than zero, ground reflections 
also contribute to the diffuse insolation on the tube array. Such re- 
flections can in some cases be a significant component of the total in- 
solation on the array, and are readily taken into account v/ith the relation 
[ 2 ] 

^ PQ(^d ^ 

where Pq is the reflectance of the ground. 

III. EFFECT OF DESIGN PARAMETERS ON COLLECTOR PERFORMANCE 
A. .Performance Index F^/F' Most of the design and operating 
characteristics of the collector are contained in the performance index 
Fg/F‘ defined by Eq. (7), Section II. A. The performance index is a function 
of two dimensionless quantities Aj and defined by Eqs. (31) and (32) 
below. 


Ai = (02/0)1 = (1 + 4 UiPi/UlP|_F')1/2 

(31) 

A 2 “ ^ 2^1 ~ ^iU|^A|^F ' /2mCp 

(32) 


The quantity Ai is termed the thermal coupling parameter. This term 
describes the effect of the heat transfer from the annulus to the delivery 
tube fluid on the performance of the collector. X 2 is called the loss/flow 
parameter since it depends on tn€ ratio of the collector loss coefficient 
to the mass flow rate of the working fluid. Notice that X 2 is also a 

- 11 - 



function of x^. 

Figure 6 is a plot of several level curves of Fj^/F' for various values 
of the parameters Xi aid X 2 . The performance index can always be improved 
(made closer to 1) if the thermal coupling parameter aj is made to approach 
a value of 1 by insulating the delivery tube. In practice this is ex- 
pensive to do and is made unnecessary by manipulating the loss/flow parameter 
X2 • 

In the performance index definition if the term is small, then 
cosh o)25i = 1 and sinh - w 2 Ci- So, 


Fp/F' - 1/(1 + ULALF72mCp) 


(33) 


Thus, as long as 1 ^ 2 ^ \ is small, (F^/F') does not depend on the thermal 
coupling term Xi no matter what its value. In principle, this approximation 
can always be made valid by increasing the mass flow rate m to offset a 
large value of the thermal coupling parameter. In practice with the present 
collector design dimensions and loss coefficient and water as the working 
fluid, the thermal coupling term Xj cannot be made larger than about 15. 

If a limit X 2 f .3 is established for the approximation to hold, then 

f 

m > 2.5 kg/hr. per collector tube. This minimum flow rate yields a tempera- 
ture rise of 15°C from 66°C at = 946 watts/m^ in one pass through the 
collector tube. Linder these operating conditions F' = .995 and Fj^ = .974. 
Higher flow rates improve Fj^ but decrease the temperature rise in one pass. 

Figure 7 shows the effect of changing the flow rate on the tempera- 
ture distribution along the delivery tube and annulus. again is 946 

watts/m^ and T^^^ = 20°C. The flow rates are 4.5 kg/hr. and 1 kg/hr. per 
collector tube. At the 4.5 kg/hr. flow rate most of the temperature rise 
in the collector takes place in the annulus. This is in agreement with 
the performance index analysis that says the thermal coupling is unimportant 
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at a sufficiently high flow rate. At 1 kg/hr., the thermal coupling is large 
with a large temperature rise in the delivery tube. As can be seen from 
Fig. 7, the collector is losing heat at a higher temperature than the 
temperature delivered at the outlet of the collector. This is not a 
particularly good operating condition even though Fp^ for this case is 
.925. The same temperature rise can be obtained by passing the fluid 
through several collector tubes in series at a higher flow rate. En- 
gineering trade-offs can be made on the basis of flow rate and tempera- 
ture rise reflected in the performance index F^^ versus pumping power and 
capacity required to push the fluid through the collector. 

B. Loss Coefficient The loss coefficient defined by Eq. (8) 

in Section II.B. is plotted in Figure 8 for two extreme ambient temperature 
conditions. The same data is presented in Tables 2 and 3 which include the 
intermediate surface temperatures T5 and Tg. Several conclusions can be 
drawn from the data in Figi 8 and Tables 2 and 3. First, since varies 
from 0.3 to 1.4 watts/m^°C in the potential operating range of the col- 
lector, the collector can be described as having a very low loss coefficient. 

, ..Secondly, decreases somewhat with decreasing ambient temperature 
because the value of U|^ is controlled by the radiation loss from the col- 
lector. This is important in heating season applications where the am- 
bient temperature is low. 

Third, increases gradually with increasing operating temperatures. 

As discussed in Section III. A, it is undesirable to operate the collector 
at a large temperature gain in one pass because of the decrease in per- 
formance index. If the collector is operated at a modest temperature gain, 

U|^ can be treated as a constant over this temperature range. Thus, linear 
theory in heat balance calculations can be used. has been treated this 
way in our computations. 
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Fourth, the heat loss through the collector is so small that the cover 
tube temperature is within a few degrees of ambient temperature even at 
high absorber temperatures. This moans that is most strongly dependent 
on the radiation loss coefficient hi. 

Last, because of the low loss coefficient, gases as well as liquids 
can be used as the working fluid. Even though the heat transfer from the 
absorber tube to the working fluid is much poorer with gases than with 
liquids, the consequent rise in absorber surface temperature does not 
increase the heat loss appreciably. Air is a particularly attractive working 
fluid because it decreases the rooftop weight of the collector and does not 
have any spill or leak problems. 

C. Insolation and Tube Spacing The components of insolation dis- 
cussed in Section II. C. can be combined into a total effective insolation 
which is used in Eq. (1) and is given by 

^eff ' ^Bh('^T V ^ i ^t) (34) 

where use has been made of Eqs. (12), (18), and (28). For simplicity, 

Eq. (34) and the discussion that follows assume that the component of 
insolation due to ground reflections is equal to S^, the uniform component 
from the sky dome; the effect of this assumption is the same as if the 
tilt s were taken equal to zero in Eq. (30), and does not materially 
affect the conclusions that follow. 

The factors Rj, W, F^y, Fj^, and F in Eq. (34) all depend on tube 
spacing d, and in each case it is found that the energy available to 
the tubes increases as d increases. In the case of Rj, the dependence on 
tube spacing enters through the shading factor g(ui) defined by Eq. (11) 
and used in approximated form in Eq. (5). The ratio Sg^/SgQ from Eq. (13) 
is shown in Fig. 9 for an equinox day and several different tube spacings 
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d. Figure 10 sliov/s this ratio for the winter and sunner solstices at 
latitude L = 40°fi and for collector tilts of 52° and 28°, favoring 
the heating and cooling seasons respectively. In Fig, 10a tube spacing 
d = 2De is assumed. Both Fig. 9 and 10 illustrate the fact that the cylin- 
drical symmetry of the tubes allows nearly all of the bean component to be 
intercepted in the early morning and late afternoon if the spacing between 
the tubes is large enough. The beam component incident on a flat surface 
at those times is generally low because of tlie cosine intercept factor, 
and consequently the operating efficiency is reduced below that for solar 
noon. The fact that an array of spaced tubes can make efficient use of the 
available light over most of a day means that the overall operating efficiency 
of a tubular array should be judged on a daily basis rather than on an 
instantaneous basis, as is frequently done with flat plate collectors. 

The factors W and F^j give rise to the dependence of the back-reflected 
beam component ^BR on tube spacing. Figure 11 shows the ratio Sg^^/pSgQ as 
a function of u for four different tube spacings d. These plots are for an 
equinox day with the plane of the screen tilted at s = L, and the screen 
distance has been taken as Dg = 3Dg/2. It is apparent that the back- 
reflected component Sg^ increases with d, and in general can be an ap- 
preciable fraction of ^30* 

Figure 12 shows the shape factor F^^ g(G) from Eqs. (24) and (29) for 
d = 2Dg and Dg = 3Dg/2, the values used with present Owens-Illinois col- 
lector arrays. Also shown in Fig. 12 is the sum F^^ g, which is seen to be 
nearly constant along the screen. This behavior is typical of other tube 
spacings of interest, and greatly facilitates the evaluation of the function 
F defined by Eq. (27). Table 4 shows several average values of F determined 
by using plots like Fig. 12, and also indicates the maximum positive and 
negative excursions from these values. The largest excursions occur for 
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small spacings, and amount to about 15% of tiie average value. For larger 
spacings, the excursions are typically under 5%. 

The effect of tube spacing on the overall thermal performance of a 
collector array is best illustrated by considering the energy outputs 
expected from arrays with different tube spacings under different operating 
conditions. Equation (34) gives the total insolation on a tube in an array 
with spacing d, while Eq. (1) gives the rate at which the tube produces 
energy. If Sp = RpSgj^ + is the insolation in the plane of the backing 
screen and d the center-to-center distance between tubes, one can define 
an array efficiency on an active installed area basis as follows; 

n(u) = Qjj/^dSp 

D4 F. 

“ (35) 

Figure 13 shows n((i>) plotted against the reduced variable (T^.^ - Tg)/Sp 
for an array in which d = 2Dg and Dg = 3De/2. For comparison, a similar 
curve for a two cover non-selective flat plate collector is 'also shown [4]. 
In arriving at these curves the following assumptions have been made. 

An equinox day has been selected at 40°N latitude, and for simplicity 
the collector tilt is taken to be s = 40°. Losses from the tubular array 
are assumed to be entirely radiative and characterized by an emittance of 
E - 0.07; the absorptance is a = 0.86, the transmittance t = 0.92, and 
the backing screen reflectance p = 0.85. These numbers represent average 
values inferred from optical measurements on the selective coating, cover 
tube glass, and backing screen respectively. In the case of the backing 
screen, the reflectance is that of an outdoor white paint on plywood. 
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In order to evaluate a day was chosen on which the total 

radiation on a horizontal surface was 20 x 10^ KJ/m^ day. The methods 
of Liu and Jordan [5] were used to divide this total into hourly beam and 
diffuse components of insolation, with the results shown in Table 5. 

The most striking feature of Fig. 13, aside from the low loss coefficient, 
is the fact that n(w) increases with hour angle w. This behavior is unlike 
that of flat plate collectors, whose efficiencies are generally independent 
of oj, and results primarily because of the back-reflected beam component 
and the symmetry of the tubes to the beam component of insolation. 

Figure 13 suggests that a better figure of merit than n(u) for 
evaluating the thermal performance of a tubular array is the daily 
efficiency n^. defined as 

eS n(w) 

% = 

(36) 

where the sums are over the hourly insolation values for a day. The 
quantity rip has been calculated for various tube spacings and operating 
temperatures, under the same assumptions used in generating Figure 13. 

T[ie results are shown in Fig. 14. It is evident that the. largest 
energy output is obtained for d = De (i.e., close packed tubes) when 
the collector is operated near ambient temperature. At higher operating 
temperatures, however, maximum efficiencies are obtained at successively 
greater spacings. Thus, at T.„ - T^ = 111°C, maximum energy output is 
obtained for d = 1.5De, while at T.„ - T, = 167°C, the maximum occurs 

in a 

at d = 2De. The basic reason for this behavior is that, on the one 
hand, the greatest ainount of radiation is intercepted, per unit of in- 
stalled area, when the spacing is small. On the other hand, the ratio 
of energy available per tube to energy lost by the tube increases with 
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spacing, so that each tube can produce more useful energy as the spacing 
is increased. 

The increase in the ratio of energy ‘ intercepted to energy lost by 
a tube with tube spacing also means that an array becomes capable of 
operating at higher temperatures as the spacing is increased. This fact 
is emphasized by Figure 15, in which the results of Figure 14 are plotted 
in a different form. This plot is similar to Figure 13, except that Hq 
instead of in((i)) is plotted on the ordinate while the reduced variable 
on the abscissa contains the daily total insolation rather than Sp. 
Although these results are based only on a single daily total insolation 
ESp = 24.5 X 10^ KJ/m?, calculations indicate that the curves can be 
applied to other total insolations with only minor error. The error, 
which is found to be only a few percent of the efficiencies shown in 
Figure 15, arises mainly from the fact that beam and diffuse insolation 
are distributed differently for different daily totals. 

An analysis of the effect of screen distance Dg on array efficiency 
can be madfe by methods which are similar to those used to analyze tube 
spacing. The results show, for example, that for d = 2De the optimum 
screen distance is Dg = SDe/Z for all operating temperatures of interest. 
More important, however, is the fact that these results indicate that 
array efficiency is insensitive to Dg over rather wide limits. In fact, 
for 0.5 5 Dg/De ^ 4, the efficiency remains constant to within 10%. 

Figure 14 shows that array efficiency also varies rather slowly with tube 
spacing. This insensitivity of efficiency to both Dg and d means that 
precise alignment of the tubes is not essential to good collector perform- 
ance, a fact which makes installation of collector arrays simpler. In 

addition, it offers array design flexibility for those situations in 
which architectural or other considerations might call for a variable 


- 18 - 



spacing and backing distance in the same collector system. 

It should be emphasized that while Figure 14 gives the optimum tube 
spacings under various conditions from the thermal performance point of 
view, economic considerations may dictate altogether different spacings 
for overall optimum system cost effectiveness. The array cost per unit 
of installed area can generally be expected to decrease with increased 
tube spacing since the hardware required decreases with increased spacing. 
If performance is viewed in terms of collector cost per unit of energy 
delivered, the optimum tube spacings are expected to be somewhat larger 
than those shown in Figure 14. 

IV. EXPERIMENTAL TESTS OF COLLECTOR ARRAYS 
Experimental tests have been made both indoors [6] and outdoors 
on a variety of tubular collector arrays using air, water, and ethylene 
glycol - waiter mixtures as the heat transfer fluid. These tests v/ere 
at ambient temperatures ranging from about -10°C to 25°C, and at operat- 
ing temperatures ranging from ambient temperature up to about 130°C. 

The tube spacing in all arrays was fixed at d = 2De, and the screen 
distance was Dg = 3De/2. Inlet and outlet temperatures were measured 
using Type T thermocouples inserted into the fluid stream at the ends of 
the manifold. For liquids, a rotameter was used to indicate flow rate, 
with more precise values being obtained by timing the accumulation of 
fluid in a graduated cylinder. An Eppley model 645-48 pyranometer mounted 
in the tilt plane of the collector was used to monitor total insolation. 
Equation (34) evaluated for d = 2Dg and Dg = 3Dg/2 would give 
for the arrays tested if the diffuse component were distributed 
uniformly over the sky dome; however, outdoor tests were usually 
made on clear days, when was probably confined to the region of the 
solar disk, and indoor tests were made with a simulator, which produces 
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no diffuse component. It should, therefore, be a valid approximation 
for all tests to write in the form 

^eff " " cosi;;)! (37) 

where is the total insolation on a horizontal surface and a is given 
by Eq. (22). Equations (37) and (35) can then be used to calculate 
instantaneous efficiencies for comparison with the test data. 

The details of the indoor simulator tests have been reported else- 
where [6]. The results are reproduced in Figure 16a, along with calculated 
curves based on the collector parameters shown in Table 6. The data 
clearly show the dependence of n(o)) on hour angle oj predicted by the col- 
lector model. In arriving at the calculated curves, it is assumed that 
the cover tube transmittance is t = 0.92; the value pA = 0.5 is inferred 
from the intercept efficiencies at the three angles, and the value U|^ = 

1.0 watt/m^°C from the slope of the data. This value of the flow 
rate (34 kg/tin. m^ of installed area), and assumed values for the heat 
transfer coefficients can then be used in Eq. (7) to determine the Fp 
shown. The absorptance a = 0.86 is then the value that best fits the 
data. 

Equation (22) and Fig. 11 indicate that a value of pA = 0.6 is to be 
expected if p = 0.85. Whether the measured value of 0.5 is lower because 
of a lower reflectance or because end effects have been neglected in the 
calculated value of a is not certain; in any event, the agreement is 
considered to be satisfactory. 

Optical and thermal testing of the tubes used in this particular test 
unit suggest an average coating emittance of e = 0.09. ihe data, on the 
other hand, suggest an effective emittance of e = 0.12 if it is assumed 
that about 10% of the experimental (J^ is due to losses in the manifold. 

The difference in these values is believed to be due to uncoated areas 


- 20 - 



near the ends of the absorber tubes and to coating scrapes along the le’ nths 
of the tubes. These scrapes occurred during the early experimental assembly 
of the test unit and have been eliminated on present tubes. Indeed, since 
the emittance of uncoated glass is about 0.9, less than 5 % of the absorber 
tube area need be scraped or uncoated to account for the above difference 
in effective emittance. 

Figure 16b shows the results of outdoor tests made at Toledo, Ohio 
(L = 42°N) on an array consisting of eight tubes. The tube axes were 
oriented north-south and the collector tilt was s = 51°. The data sum- 
marize tests made during the period January 23, 1975 to February 14, 1975 
when the declination 6 varied from -20°C to -14°C and tlie ambient temperature 
from -10°C to +10°C. Wind speed was not recorded. Pure ethylene glycol 
was used as the heat transfer fluid at all temperatures (up to 130°C). The 
average flow rate was about 25 kg/hr. m^. Also shown are calculated curves 
for w = 0° and oj = 40°, based on the parameters shown in Table 6. In this 
case the angular dependence of n(w) is not as sharply defined by the data 
as in Figure 16a, but it is reasonable to suppose that this fact is due 
largely to the variability of outdoor conditions, which are not as well 
controlled as indoors. In any event, the calculated curves are consistent 
with the data. For these tests the parameter pA = 0.73 is reasonably 
consistent with the expected value of pA = 0.6. As in the case of the 
indoor tests, the difference could be due either to screen reflectance or 
the neglect of end effects in the calculated value. The loss coefficient 
U|^ is somewhat lower than for the unit tested indoors, the improvement prob- 
ably being the result of the elimination of coating scrapes on the tubes 
used in the present unit. If lOi^ of the experimental Uj^ is assigned to mani- 
fold losses, it is found that the effective emittance is e = 0.09, as com- 
pared with an expected emittance of e = 0.07 based on optical tests of the 

- 21 - 



coating. The difference is probably caused in part by the fact that, even 
though the tubes had no coating scrapes, the ends of the absorber tubes were 
still uncoated and therefore, had the emittance of bare glass. Since recent 
improvements in coating procedure now allow the entire absorber tube to 
be covered, future collector arrays are expected to have even lower loss 
coefficients. 


V. TRANSIENT BEHAVIOR 

The general time and space dependent solution to Eqs. (2) and (3) 
is still under investigation. However, a solution has been obtained for 
the case where thermal coupling between the two fluid passages is un- 
important. This situation corresponds to the operating conditions already 
discussed in Section III. as being desirable for good performance. 

If the thermal coupling in Eqs. (2) and (3) is neglected anc! the 
boundary conditions are assumed to be 

Ti(0,t) = T.^(t) and Ti(?i,t) = T2(5i,t) (38) 

it can be shown that for t > (a + l)ci, when start-up effects are no 
longer important, 

% = mCp{To(0,t) - T.^[t-(a+l)?i]} 


= Fp(l)A^ {ax<Sgff> - \ ^ [T.^[t-{a+l)ei]-Tg]} 


(39) 


Here, Fr( 1) is given by Eq. (7) with the coupling parameter Ai equal to one 
and <Sg^^> is defined by 


F- p 

[TTiTj 

Jo 


^^eff>" Fp(l)ci I Sg^^(t-0)e 


(40 
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The difference between the times at which the inlet and outlet temperatures 
are evaluated is equal to the time required for an element of fluid to 
move through the collector. The quantity is a weighed average of 

the effective insolation over the time that the fluid is in the annulus. 

The time delays are of interest in the tubular collector because under 
normal or "good" operating conditions, an element of water can be resident 
in the collector for about thirty minutes, and this situation has im- 
portant implications for the control characteristics of the collector. 
Equation (39) is an exact solution for a well-insulated delivery tube 
or a single pass collector, and should be a good approximation for the 
Owens-Illinois tubular collector as long as thermal coupling effects are 
not significant. In fact, preliminary tests during which T^.^ was allowed to 
vary with time indicate that Eq. (39) properly accounts for time (indicate 
that Eq. (39) properly accounts for time) delay effects. 

VI. DISCUSSION 

TM 

The Owens-Illinois, Inc. SunPak ^ solar collector described in 
this paper incorporates several desirable collector features. Collector 
evacuation and selective coating produce an extremely low heat loss coef- 
ficient Uj^. Some aspects of the low loss coefficient are discussed in 
Section III.B. Further, the low heat loss enables the collector to 
operate with good efficiency at temperatures high enough to drive existing 
heating and cooling devices presently using fossil fuel energy sources. 

This means that retro-fitting solar energy to existing installation is 
now feasible. In addition, present design heating and cooling devices 
can be used without de-rating their capacities due to inadequate 
temperatures obtained from the collector. Alternatively, and possibly 
more important, the SunPak collector can operate with good efficiency 
at more moderate temperatures on days when the insolation is low. In fact. 
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for a collector temperature of 65°C and an ambient temperature of 0°C, only 
80 v/atts/m^ of insol atic;) are required in order to obtain useful energy 
from the collector. 

The performance index (Fj^/F') described In Se'^tion III. A. can be 
made greater than .9 for all applications envisioned to date. The thermal 
coupling between the dt^li very tube and the annulus does adversely affect 
the performance index if the operating conditions are not carefully 
examined. However, suitable flow rates and array configurations can be 
designed so that the thermal coupling becomes negligible. 

As discussed in Section III.C., the effective insolation on a tube, 
has a time dependence that differs from that of most collectors. 
Because of the tube-to-tube spacing, the collector intercepts direct 
beam radiation uniformly for eight hours with the ,^resent collector 
design. In addition, the back reflective screen pe)mits the recovery of 
about fifty per cent of the radiation that falls between the tubes. This 
means the instantaneous efficiency n(w) as usually defined increases 
toward the beginning and end of the day for a south facing array in the 
northern hemisphere. As a consequence, the collector array performance 
cannot be realistically evaluated on the basis of instantaneous efficiency, 
but rather must be judged on the basis of efficiency or energy output over 
periods no shorter than a day. 

The insensitivity of collector efficiency to operating temperature, 
ambient temperature, and wind speed has already been mentioned. An impor- 
tan consequence is flexibility in system design, since temperature 
conditions can vary over rather wide limits with no serious degradation 
of collector performance. This insensitivity also means that the collector 
can operate with air as the heat transfer fluid, with nearly the same 
efficiency as if a liquid were used. Because of its low film coefficients. 
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air usually imposes a severe penalty on flat plate collectors by causing 
the absorber surface to run hotter than with a liquid. In the case of 
the evacuated tubular collector, the absorber surface runs hotter than 
with a liquid, but Figure 16 shows that no significant decrease in ef- 
ficiency will result. No modification of the basic tubular collector 
elements is necessary for operation with air. All that is necessary, 
in fact, is a larger delivery tube and a different manifolding system 
designed to give the proper air pressure drops through the collector 
array. Preliminary testing of air operated tubular arrays at the 
Owens-Illinois test site confirms that the performance is essentially 
as good as with a liquid heat transfer fluid. 

Although not all of the predictions of the collector model developed 
in this paper have been verified by direct test data (for example, arrays 
with tube spacings other than d = 2De have not been tested,) it is believed 
that the data presented in Section IV. substantiate the most important 
features of the model. The tube spacing d = ZDg used with present Ov/ens- 
Illinois collectors was selected on the basis of both thermal performance 
and cost-effectiveness of the collector. Figure 15 shows that for this 
spacing the collector is less efficient at low temperatures (or high in- 
solations) than if the tubes were more closely spaced. While the tubes 
could be more closely spaced to give better thermal performance at low 
temperatures, performance at high temperatures (or low insolations) would 
suffer and it is not clear that the collector would be as^ cost effective. 

Finally, the collector tube material, glass, is already known to have 
good weathering and chemical durability properties, assuring long life of 
the collector. The present production collector is sold as a unit con- 
sisting of the collector tubes and manifold which has headers incorporated 
in it. The units can be installed by one man either on the rooftop or 
other collection site. No mechanical lifting equipment is required for 
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installation. Large arrays can be built up from the individual units 
with little difficulty. 
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Delivery Tube 

I.D. 

- Di 

9 

mm 

Delivery Tube 

O.D. 

- D2 

- 12 

mm 

Absorber Tube 

I.D. 

- O 3 

- 39 

mm 

Absorber Tube 

O.D. 

- D 4 

- 43 

mm 

Cover Tube 

I.D. 

- Ds 

- 49 

mm 

Cover Tube 

O.D. 

- Ds 

- 53 

mm 

Active Tube Length 

- £ 

- 1067 

mm 

Subscripted parameter val 

lues ' 

in the body of the 


paper refer to the numbering scheme established 
here. 


Table 1: Collector tube dimensions and subscript 

identification 
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Absorber Temperature 


Inside & Outside Cover Temi 


Loss Coefficient 
Watts 

0.286 
0.303 
0.321 
0.340 
0.360 
0.381 
0.402 
0.425 
0.449 
0.473 
0.499 
0.526 
0.554 
0.583 
0.613 
0.645 
0.677 
0.711 
0.746 
0.782 
0.820 
0.859 
0.899 
0.941 
0.983 
1.028 
1.074 
. 1.121 
1.170 
1.220 


Table 


®C 

0 . 

10.0 

20.0 

30.0 

40.0 

50.0 

60.0 

70.0 

80.0 
90.0 

100.0 

110.0 

120.0 

130.0 

140.0 

150.0 

160.0 

170.0 

180.0 

190.0 

200.0 
210.0 
220.0 

230.0 

240.0 

250.0 

260.0 

270.0 

280.0 
290.0 


2: Lbss coefficient for 


“C 

«c 

- 19.9 

- 19,9 . 

- 19.8 

- 19,8 . 

- 19.7 

- 19,7 

- 19.6 

- 19,6 

- 19,5 

- 19.5 

- 19.4 

- 19.4 

- 19.3 

- 19.3 

- 19.1 

- 19,2 

- 19,0 

- 19.1 

▼ 18.8 

- 18.9 

- 18.6 

- 18.8 

- 18.4 

- 18.6 

- 18.2 

- 18.4 

- 18.0 

- 18.2 

- 17.8 

- 18.0 

- 17.5 

- 17.7 

- 17 . 2 . 

- 17.5 

- 16.3 

- 17.2 

- 16.6 

- 16.9 

- 16.3 

- 16.6 

- 15.9 

- 16.3 

- 15.5 

- 15.9 

- 15.1 

- 15.6 

- 14.7 

- 15.2 

- 14.2 

- 14.7 

- 13.7 

- 14,3 

- 13.2 

- 13.8 

- 12.6 

- 13.3 

- 12.1 

- 12.8 

- 11.5 

- 12.2 


, = -20°C 

a 
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Loss Coefficient 
Watts 
M^~^C 

Absorber Temperature 
°C 

Inside & Outside Cover Temp. 
°C °C 

0.503 

50.0 

40.1 

40.1 

0.527 

60.0 

40.2 

40.2 

0.553 

70.0 

40.4 

40.3 

0.579 

. 80.0 

40.5 

40.4 

0.607 

90.0 

40.6 

40.6 

0.636 

100.0 

40.8 

40.7 

0.665 

110.0 

41.0 

40.9 

0.696 

120.0 

41.2 

41.1 

0.729 

130.0 

41.4 

41.3 

0..762 

140.0 

41.6 

41.5 

0.796 

150.0 

41.9 

41.7 

0.832 

160.0 

42.1 

41.9 

0.869 

170.0 

42.4 

42.2 

0.907 

180.0 

42.7 

42.4 

0.947 

190.0 

43.0 

42.7 

0.988 

200.0 

43.3 

43.0 

1.030 

210.0 

43.7 

43.3 

1.074 

220.0 

44.1 

43.7 

1.119 

230.0 

44.5 

44.0 

1.165 

240.0 

. 44.9 .. 

44.4 

1.213 

250.0 

45.4 

44.8 

1.263 

260.0 

45.9 

45.3 

1.313 

270.0 

46.4 

45.7 

1.366 

280.0 

46.9 

46.2 

1.420 

290.0 

47.5 

46.7 


Table 3: Loss coefficient Uj^ 

for T, = 40°C 

a 
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Spaci^ /> d 


De 


0 


l.BDe 

0.14 

± 

0.02 

2De 

0.343 

+ 

0.013 

3Dg 

0.799 

+ 

0.014 

4Dg 

1.280 

+ 

0.014 


Table 4: Values of F [Eq. (27)] for different tube spacings d. 
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Hour From Solar Hoon 


Beam Component 

Oiffuse Component 


^BH 

^BP 

^BO 

Sd 


(W/m2) 

(W/m^i 

(W/m2) 

(W/m2) 

± 1/2 

606 

792 

798 

174 

i 1-1/2 

543 

707 

767 

158 

± 2-1/2 

442 

577 

726 

139 

± 3-1/2 

315 

41 0 

675 

104 

± 4-1/2 

155 

202 

527 

66 

± 5-1/2 

38 

50 

382 

22 

^d^ ~ 

19.9 X 1Q3 

kJ/m^ day 




24.5 X 1Q3 

kJ/m^ day 



Equinox day (6 

= 0) at L 

= 40°N, s 

II 

o 

o 



Table 5: Hourly beam and diffuse insolation on 

equinox day 


BB££eS>im £AC£i IB 


OF tOOR QUALOT 
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Test 


F,^ a , T pA (W/m2“C) 


Indoors 


(NASA LeRC) 

0.975 

0.86 

0.92 

0.5 

1.0 

Outdoors 

(0v.’ens-Illino1s) 

0.98 

0.85 

0.92 

0.73 

0.85 


Table 6; Collector parameters derived from test data 
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Figure 1: 
Figure 2: 
Figure 3: 
Figure 4: 
Figure 5: 
Figure 6: 

Figure 7: 

Figure 8: 

Figure 9; 
Figure 10: 

Figure 11; 
Figure 12: 

Figure 13: 

Figure 14: 
Figure 15: 
Figure 16: 


Figure Captions 

Details of collector tube assembly. 

Collector array at Owens-Illinois test site. 

Control volumes for heat balance of Eqs. (2) and (3). 

Thermal loss network for collector tube assembly. 

Geometry of a spaced tube collector array. 

Performance index level curves for different coupling and 
loss-flow parameters. 

Temperature rise along a single collector tube for two dif- 
ferent fluid flow rates. 

Dependence of the loss coefficient on absorber tube surface 
temperature and ambient temperature. 

Dependence of Sgg on time for different tube spacings. 

.Dependence of Sgp on time at the winter and summer solstices, 
for two different collector tilts. A spacing d = 2De is assumed. 
Dependence of Sgg on time for different tube spacings. 

Shape factors F^^ ^(6) for d = 2De and Dg = 3D6/2- The sum 
^dx s nearly constant across the screen. 

Instantaneous efficiency of a two cover, nonselective flat 
plate collector and a south facing tube array with d = 2De. 

An equinox day and s = L are assumed. 

Dependence of daily efficiency of a tubular array on tube spacing 
for different operating conditions. 

Dependence of daily efficiency of a tubular array on (T.^-Tg)/ESp 
for different tube spacings. 

Measured and calculated efficiencies of two tubular test arrays, 
(a) Data from NASA LeRC indoor simulator, (b) Outdoor data 
taken at Owens-Illinois test site. 
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NOMENCLATURE 


The following definitions of parameters and variables are used through 
out the paper. We have attempted to use terms that are in current use 
in the body of literature on solar energy where possible. Subscripts 
are used in the paper to denote which component the parameter applies. 

A “ Surface area 
A^ - Cross-sectional area 
Cp - Specific heat 
D - Diameter 

Dg - Distance from tube axis to reflecting screen 

E - Black body emissive power 

F - Geometric shape factor 

Fp - Collector performance index 

F' - Collector plate factor 

J - Radios! ty 

L - Latitude angle 

P - Perimeter 

Q - Heat rate 

S - Insolation 

T - Temperature 

U - Overall heat transfer coefficient 
d - Center line to center line tube spacing 
h - Convective heat transfer coefficient 
k - Thermal conductivity 

£ - Collector tube length 

• 

m - Mass flow rate 
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s - Collector tilt angle to horizontal 
t - Time coordinate 

A 

t - Reduced time coordinate 

X - Space coordinate 

K - Reduced space coordinate 

a - Collector absorptance 

e - Emittance 

i5 - Declination angle 

A - Back-reflected light parameter 

p - Reflectance 

Pf - Heat transfer fluid density 
w - Hour angle 
T - Collector transmissivity 
Xi - Thermal coupling parameter 
X 2 - Loss/ flow parameter 
a - Stefan-Boltzmann constant 
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